Kinetic Magnetism and Orbital Order in Iron Telluride 



Ari M. Turner, 1 Fa Wang, 1,2 and Ashvin Vishwanath 1,2 

1 Department of Physics, University of California at Berkeley, Berkeley, CA 94720 
Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720 

(Dated: May 23, 2009) 

Iron telluride (FeTe), a relative of the iron based high temperature superconductors, displays 
unusual magnetic order and structural transitions. Here we explore the idea that strong correlations 
may play an important role in these materials. We argue that the unusual orders observed in FeTe 
can be understood from a picture of correlated local moments with orbital degeneracy, coupled to 
a small density of itinerant electrons. A component of the structural transition is attributed to 
orbital, rather than magnetic ordering, introducing a strongly anisotropic character to the system 
along the diagonal directions of the iron lattice. Double exchange interactions couple the diagonal 
chains leading to the observed ordering wavevector. The incommensurate order in samples with 
excess iron arises from electron doping in this scenario. The strong anisotropy of physical properties 
in the ordered phase should be detectable by transport in single domains. Predictions for ARPES, 
inelastic neutron scattering and hole/electron doping studies are also made. 



PACS numbers: 

I. INTRODUCTION 

The discovery of high temperature superconductivity 
in a class of iron based materials^ has opened a new 
route to high temperature superconductivity besides the 
ones operating in the copper oxide materials. Follow- 
ing the initial discovery in LaOFeAs (1111 materials), a 
number of new classes of materials were discovered that 
shared similar properties, notably the 122 materials (like 
BaFe2As2). In these systems, (collectively referred to as 
the FeAs materials) the undoped compound is a metallic 
SDW system, with ordering wavevector (tt, 0), which on 
doping leads to a superconducting state. The magnetism 
in these materials is believed to arise from Fermi sur- 
face nesting, given the presence of an electron and hole 
pocket separated by (tt, 0) in the LDA band structure 
calculations of these materials. In fact, theoretical stud- 
ies had predicted this ordering before it was confirmed 
in neutron scattering experiments. Moreover, the order- 
ing moment is typically quite small, eg. of order O.SfiB in 
LaOFeAs, and the absence of a Curie- Weiss form of mag- 
netic susceptibility above the ordering temperature have 
been invoked as evidence for the itinerant character of 
the magnetism. Finally, signatures of an excitation gap 
appear in optical conductivity experiments, on cooling 
through the SDW transition^ 

An important recent development has been the dis- 
covery of a new class of materials FeSe and FeTe, which 
share the square lattice Fe structure of the FeAs ma- 
terials and are believed to be closely related. Indeed, 
superconductivity has been observed in FeSe even in the 
absence of doping at 8K, rising to 37 K on application 
of hydrostatic pressure. The chemical simplicity of these 
materials, as well as the absence of a pnictide group ele- 
ment, may offer valuable clues to isolating the physics of 
the iron based high temperature superconductors. One 
notable difference from the FeAs materials though, is in 
the magnetism. While FeSe is non-magnetic even at sto- 



ichiomctry, the FeTe materials are magnetically ordered 
metals, but with a more complicated kind of order than 
seen in FeAs, shown in Fig. [TJ The ordering wavevec- 
tor is (7r/2, tt/2), in contrast to the (tt, 0) ordering of the 
FeAs compounds and also the (tt, tt) ordering of the insu- 
lating parent compound of the cuprate superconductors. 
(Note the wavevectors here are defined with respect to 
an unfolded zone comprising of a unit cell of one iron 
atom, oriented along the iron square lattice. The actual 
basis vectors are x ± y because of the alternating po- 
sitions of the tellurium ions; therefore, crystallographic 
studies use a doubled-unit cell, and a and b are along 
the diagonals of the iron lattice.) The order sets in via 
a first order transition at 87K , and is accompanied by 
a monoclinic distortion-. In the presence of excess iron, 
i.e. Fei+yTe, the commensurate order is found to evolve 
into an incommensurate spiral. Whether this complex 
magnetic behavior is important to understanding other 
physics in this material is presently unclear. So far, it is 
unique to FeTe, where superconductivity only occurs on 
substantial alloying with sulfur FeSo^Teo.s or selenium 
FeSeo.$Teo.5. However, besides the interest in under- 
standing the origin of this unusual magnetism, there are a 
number of indicators that point to the presence of strong 
correlations in FeTe, which would be an important fact 
to establish in these materials. (1) The ordered mag- 
netic moment in FeTe is large, ~ consistent with 
a localized S = 1 at every site. (2) Above the order- 
ing temperature the magnetic susceptibility falls off in a 
Curie Weiss fashion, roughly consistent with the observed 
ordered moment and transition temperature. (3) given 
the absence of Fermi surface nesting at this wavevector, 
a spin density wave scenario seems less favorable than 
a local moment picture. Furthermore, ARPES experi- 
ments on this material see no obvious nesting at the de- 
sired wave-vector—, and (4) optical conductivity, which 
observes a clear SDW gap in the FeAs materials, does 
not see an analogous gap in FeTe 5 . Another mechanism 
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FIG. 1: The commensurate (■§,■§) spin ordering pattern of 
Fei+jjTei, close to x = 0. The spins reside on the iron atoms, 
which form a square lattice. The tellurium atoms alternate 
above (x's) and below (o's) the iron plane. When the spins or- 
der, the square lattice of irons distorts into an approximately 
rhombic lattice. We assume that the monoclinic distortion 
is a consequence of the spin ordering. The tellurium atoms 
move towards (or away from) the exceptional spin in each 
parallelogram. Thus, the planes of tellurium atoms shift in 
opposite directions, producing the monoclinic distortion. 

is therefore needed to explain the peculiar spin ordering. 
Here we will assume that the magnetism in FeTe arises 
from local magnetic moments on the iron sites, that arise 
from strong correlations. However, given that FeTe is 
metallic (with a residual resistivity of 0.2mf2 cm) - we 
will have to consider them as being 'self-doped'. The 
key point of this paper is that several of the puzzling 
magnetic properties of FeTe can be naturally explained 
if we assume it is near a correlated insulating state with 
spin and orbital degeneracies. Furthermore, the struc- 
tural distortion here will be explained as arising, at least 
partially, from orbital ordering - rather than spin lattice 
interaction as is usually assumed. A model for FeTe's 
ordering in which the lattice distortion is assumed from 
the outset is in e . Other scenarios based on an itinerant 
electron picture have also been put forward^, such as the 
possibility that electron doping is large enough to change 
the Fermi surface shape and lead to nesting 8 . 

The main assumptions we make in this paper are: 
(i) FeTe is proximate to an S — 1 magnetic insula- 
tor. (ii)Each site of this insulator has orbital degener- 
acy (g?i, (I2) and the Jahn Teller effect leads to orbital 



ordering which orients the orbitals along the diagonal di- 
rection of the Fe lattice. The orbital e?i accommodating 
the local moment is uniformly oriented on all sites, (in) 
A small excess density of charge carriers is present in the 
other Jahn Teller orbital d^- 

From these assumptions we show that the (tt/2,tt/2) 
magnetic order can be naturally explained, as well as 
the incommensuration induced by excess iron. The key 
mechanisms are the formation of one dimensional chains 
induced by the Jahn- Teller ordering, which are coupled 
together by double exchange. Part of the structural tran- 
sition in this scenario is caused by orbital rather than 
spin ordering. Since this proposal invokes an intertwining 
of the spin, charge and orbital physics, several testable 
consequences also emerge for transport and ARPES ex- 
periments. For example, transport within a single do- 
main is expected to be highly anisotropic, with larger 
conductivity along the jerromagnetically ordered diago- 
nal directions. This can potentially be probed by optical 
conductivity experiments. Indeed, a number of anoma- 
lies are already observed in transport properties, albeit 
in multidomain samples. The electrical conductivity rises 
sharply below the ordering transition while the Hall effect 
abruptly changes sign. 

Recently, a similar scenario has appeared to explain the 
magnetism and structural transition in FeAs^ ' 10 ! 11 . In 
contrast to our assumption (ii), the Jahn- Teller orbitals 
are assumed to order along the lattice directions. The 
spin wave spectrum experimentally observed in FeAs was 
argued to be much better described by such a models. 
Here we focus on FeTe, for the reasons described previ- 
ously, but our approach is very similar in spirit. If indeed 
this mechanism is more general, it allows us to unify phe- 
nomena across this family of compounds. 

This paper is organized as follows. We first describe 
microscopic strongly correlated models for FeTe, which 
have S=l and orbital degeneracy. We then consider how 
orbital ordering of a particular kind can drive magnetism, 
leading to the observed spin order. The characteristics of 
spin wave dispersion within this scenario are presented, 
and the origin of incommensurate magnetism with excess 
iron, is discussed. Finally, experimental consequences of 
this scenario for nonmagnetic properties, like conductiv- 
ity and ARPES, are described. 



II. FeTe: A STRONGLY CORRELATED 
VIEWPOINT 

A. Microscopic Model 

We first model FeTe in terms of a nearby correlated in- 
sulating state. The modifications required to account for 
metallic conduction are discussed later. We demand that 
the insulator carries net spin 5=1, per Fe atom, based 
on the ordered moment observed at low temperatures. 
Furthermore, we will require that they exhibit orbital de- 
generacy. The only pair of d orbitals that are degenerate 



FIG. 2: Two scenarios for the level spacing of the iron atoms' 
d-orbitals which lead to orbital degeneracy and S = 1. Each 
set of orbitals grouped together is assumed to fill up accord- 
ing to Hund's rule before any electrons are added to the next 
group of orbitals. This follows from assuming that the crystal 
field splitting between orbitals in different groups is greater 
than the Hund's coupling. Scenario 1 results in a single elec- 
tron having to decide between two degenerate orbitals while 
Scenario 2 results in a single hole which is orbitally degener- 
ate. We mainly discuss Scenario 1. 



in this tetragonal structure are the d xz and d yz orbitals. 
Hence we will require an odd filling of this orbital pair. 
Two possible microscopic scenarios for the d 6 configura- 
tion of the Fe +2 ion are sketched in Figured] In the first, 
there is one electron available to occupy the two degen- 
erate orbitals, while the orbital d xy is singly occupied. 
Interestingly, this is the orbital configuration suggested 
by the crystal field splitting of the Fe sites. For a perfect 
tetrahedral arrangement of Te ions, the d x 2_ y 2 and d z 2 
orbitals lie below the triplet of d xyi d xz , d yz . The distor- 
tion of the tetrahedron in this material brings the d xy be- 
low the degenerate d xz , d yz - leading to the orbital struc- 
ture in Figured Note, the sense of the distortion in FeTe 
is opposite to that in FeAs, where a similar exercise leads 
to a different orbital ordering 3 -. Since such a local picture 
of electronic orbitals may not capture the physics of FeTe 
we also point out a different scenario (Figure dt 1 ) , where 
the orbital ordering is closer to what is predicted by the 
band structure calculations in FeAs^. If we order the 
dispersing bands by their center of mass, we end up with 
the ordering shown. Here too, the electron assignment 
can lead to a S = 1 orbitally degenerate configuration, 
but now the degenerate pair of orbitals contains three 
electrons. We note that in both scenarios, the active or- 
bitals are d xyi d XZl d yzi which are also the ones expected 
to be present at the Fermi energy from weakly correlated 
band structure calculations of these materials 1 - 113 . The 
two scenarios are particle hole conjugates of one another, 
if one focuses on the active triplet of orbitals. Henceforth 
we will assume Scenario 1 for concreteness, as it corre- 
sponds better with experimental facts. The results there 
can be easily transcribed to Scenario 2, by a particle hole 
transformation. 



FIG. 3: Degeneracy lifting by Jahn- Teller distortion. The 
tellurium atoms projected into the xy-pla,ne form a square in 
the undistorted compound, which distorts into a rhombus, so 
that A > B, lowering the energy of the shaded orbital. 



B. Orbital and Magnetic Order 

Let us discuss first the the pair of degenerate orbitals 
in Scenario 1, filled by a single electron. Later we will 
include the third orbital. It is convenient to rotate or- 
bitals by 45°, and define another basis that point along 
the diagonals dx z /Yz = -j^{d xz ±d yz ). Let e?{, d\ be the 
creation operators for electrons in these diagonal orbitals. 
The Hamiltonian for this system is: 



H = H 



KE 



H„ + H. 



JT 



(1) 



where the first term is the hopping hamiltonian, and the 
second and third terms refer to interactions and coupling 
to lattice phonons that lead to the Jahn- Teller effect. 



H, 



^U(n r - l)n r - JnSir ■ S 2r 



(2) 



where n r is the electron density at site r and S a r = 
\d\ r ad ar is the spin on site r in orbital a. In the limit of 
strong repulsion U and a single electron per site, (n) = 1 
we obtain an insulating state with orbital degeneracy. 
This degeneracy is typically resolved by the Jahn- Teller 
effect. A lattice distortion which breaks symmetry and 
splits the degeneracy occurs. The precise distortion that 
is realized is hard to predict, so we will assume that it 
is indeed of the type required to obtain the structural 
transition seen in this material. This involves a uniform 
orthorhombic distortion that changes the relative lengths 
of the two diagonal bonds. If we denote by A and B, the 
classical bond lengths for the two diagonals (so A = B in 
the tetragonal state), then we will assume that the lattice 
coupling is given by: 



H 



JT 



Y{A-B)J2(d\ r d lr 



-d\ r d 2r ) + l-{A-Bf. (3) 



(See Fig. [3J) The orbital ordering implies that, in sce- 
nario 1, the single electron on each site always occupies 
the same orbital, i.e., ri\ r = 1 and n 2r — 0, or vice versa. 
For concreteness let us suppose that the X diagonal ex- 
pands. Then the electrons occupy orbital 1 (dx z ) while 
orbital 2 (dy z )is empty, because it is higher in energy. 
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Since each site has an unpaired electron, we can now de- 
rive the Hamiltonian governing their magnetic moments. 

In the insulating limit, the magnetic interaction is gen- 
erated by virtual hopping of electrons. To proceed, we 
need to specify Hke- Clearly, given the geometry of the 
Xz, /, Yz orbitals, hopping along the diagonals will be 
very anisotropic. We denote the X (Y) diagonal hop- 
ping of the dxz (d>Yz) orbital by t2, and let t' 2 be the 
hopping of each in the orthogonal direction (see Fig. rj|. 
The figure suggests that t% S> t' 2 (seeM), and for FeAs 
the nearest neighbor hoppings are comparatively small 
as well. Hence, we will simply work with the exchange 
interaction induced by t%. Note, the next neighbor hop- 
ping only operates within a single sublattice (labeled A 
and B in the figure), so here we consider just the A sub- 
lattice. The t2 hopping will introduce antifcrromagnetic 
exchange, but only along the X diagonal Jix = t^/U. 
This will lead to antiferromagnetic order along this di- 
agonal direction. Note, assuming Scenario 1, this an- 
tiferromagnetic direction will be the expanded diagonal 
of the distorted compound, which is consistent with the 
observed wavenumbers of the magnetic and structural 
distortions^. Note, although the diagonal chains are or- 
dered, there is negligible coupling between chains at this 
moment (because of the smallness of t' 2 ). Below, we will 
see there is a more important mechanism that can lock 
the magnetic order in the chains together - the double 
exchange interaction. 



C. Double Exchange 

Hopping of electrons between the antiferromagnetic 
chains can readily occur if they occupy the dy z orbital. 
However, in a orbitally ordered insulator, these are as- 
sumed to be empty. Given the empirical fact that FeTe 
is a metal, we assume a small occupation S of electrons in 
this orbital. This can arise because of 'self-doping', i.e. 
if the orbital dispersion cannot be completely neglected, 
a fraction of carriers from one of the 'filled' bands (for 
example d xy ), could be transferred to an empty band, i.e. 
this orbital. This also leads to a conducting state at sto- 
ichiomctry, as demanded by experimental observations. 
Furthermore, FeTe always occurs with a slight excess of 
Fe i.e. Fei+^Te. The electrons from the additional Fe +2 
ions also contribute to S. 

Such an excess carrier density will lead to ferromag- 
netic interactions between chains, via the double ex- 
change mechanism. An electron hopping in this nearly 
empty orbital will be Hund's coupled to the dx z elec- 
tron according to Equation [2] The Hund's coupling is 
typically large and will force both electrons to have the 
same spin. If the electron is to hop to the neighboring 
chain, the spins must be parallel. Then, it can enjoy a 
lowering of kinetic energy by —2t2- Thus, a ferromag- 
netic arrangement of spins will have a lower energy than 
an antiferromagnetic arrangement by an amount 2|t2|<5, 
which can roughly be viewed as a ferromagnetic coupling 



along the Y diagonal (strictly speaking this is a nonlocal 
interaction, and cannot be assigned solely to the diago- 
nal bond). Thus J^y — —1\ti\6. With this combination 
of exchange constants, the magnetic ordering on a sin- 
gle sublattice is shown in Figure [4] Note that it has the 
wavevector (7r/2,7r/2), as required. 

D. Coupling the Sublattices 

So far, the two sublattices (^4 and B) are independent. 
When the doping is small, the dominant spin interaction 
between nearest neighbor sites will arise from antiferro- 
magnetic exchange from the electrons in the half filled 
d xy orbital. (The nearest neighbor hopping of electrons 
in the dx z orbitals would induce a smaller ferromagnetic 
interaction controlled by the Hund's coupling, since the 
dominant hopping is expected to switch the orbital in- 
dex.) We will also invoke a coupling to the lattice to 
generate a biquadratic interaction term in order to lock 
in the commensurate wave number, leading to the net 
Hamiltonian describing the interaction of spin 1 atoms 
on a square lattice: 

H = J2 JijSi ■ Si - ^^(S, ' S j) 2 ( 4 ) 

ij (ij) 

J2X = t%/U, J2Y — — 2|i 2 |<5 and J\ > 0, and J\ is the 
same for both nearest neighbor bonds. The phase dia- 
gram of this model as a function of increasing K is in- 
cluded on the y-axis of Fig. O The Sj • Sj terms alone 
would lead to an incommensurate spiral stated, where 
the nearest neighbor spins in the spiral state are close to 
being orthogonal to one another (the "single spiral state" 
below). This state, even in the small incommensuration 
limit, is significantly different from the experimentally 
observed collinear state. 

To stabilize the commensurate , — 5) state, we take 
into account the biquadratic spin interaction (— K(Si ■ 
S2) 2 ) which is a well-known consequence of spin-lattice 
coupling. This term prefers collinear magnetism. As we 
will see below, even modest values of the spin phonon 
coupling K can induce locking of the commensurate, 
collinear phase observed in experiments. In particular, 
we show below that the critical coupling required to in- 
duce collinear order K > J 2 / J 2 x can be parameterically 
smaller than Jj . Coupling of spin and lattice is presum- 
ably essential to getting a commensurate state at this 
wavevector. 

Phase Diagram: The phase diagram can be obtained 
by taking the ansatz of a pair of coplanar spirals on the 
two sublattices with an arbitrary phase <fi between them, 
and wavevector k along diagonal X. This is probably suf- 
ficiently general to capture the ground states of Eqn. |H 
because the ferromagnetic coupling along the Indirection 
prevents the spin from varying in that direction. Let the 
azimuthal angle of spin i be 9i. Then the ansatz reads: 

i = k{x i -y i ) + (-l)*< + «£ (5) 



a) 




b) 




FIG. 4: How the spin-order on sublattices arises out of the orbital order. The solid and dotted lines show the projections of 
the dxz and dy z orbitals. a) The second nearest neighbor hopping amplitudes, which determine the order in the sublattices. 
Electrons stay in the same orbital when they hop to the second nearest neighbor. The hopping for each orbital is anisotropic, 
with an anisotropy that depends on the orbital, so each orbital forms a set of one-dimensional chains. Electrons can hop 
between nearest neighbors as well (not shown), b) Generating the spin order in the A-sublattice. (The iron atoms in the 
B-sublattice are indicated by x's.) The dxz are shaded to indicate that they are filled and the dyz orbitals are lightly doped. 
The dxz orbital form one-dimensional Mott insulators parallel to the X-axis with antiferromagnetic coupling Jix- The dyz 
orbitals form one-dimensional metals parallel to the Y direction. Ferromagnetic order along Y lowers the kinetic energy of the 
metals by about t2&. 



Here,(fc, —k) is the wave number of the spin arrangement. 
(This state is ICA, from Ref.— .) Note, the experimentally 
observed collinear state corresponds to k = tt/2 and <j) = 
tt/2. Along any row, the angle between adjacent spins 
alternates between k — </> and k + cj>. For k = <j) = n/2, 
the spins therefore alternate from parallel to antiparallcl. 
Minimizing the energy per site, 

E[k, 4>] — J 2 x cos 2k + 2J\ cos k cos <j> 

-K(l + cos 2k cos 20) - 2\t 2 \6 

(6) 

over cj) and k shows that there is a transition between a 
pure spiral with (f> — and incommensurate wavevector 
k = § + 0(t~") an d the collinear state k = tt/2, <fi = tt/2 
when AK(J 2X — K) = J\. For small K < J 2X , the 
critical K is ps J\I^J 2 x- 

We can understand this result intuitively as follows. 
The spiral order occurs to lower the nearest neighbor 
energy. The coupling of a spin to its neighbors to the left 
and right cancels in the perfectly collinear state. To take 
advantage of Ji, the neighbors should therefore make an 
angle less than 180°, and the central spin should point 
opposite to the sum of the two neighboring spins. If J\ 
is small compared to J 2 x , then the neighbors are nearly 
antiparallcl, so J\ has a very weak effect, explaining why 
the critical value for K is not of order Jj as one might 
have expected, but rather second order in Jy, 

Spin Waves: To contrast this scenario with others that 
predict the same magnetic ordering pattern, we calcu- 



late the spin wave spectrum for the model [4] We ex- 
pect inelastic neutron scattering experiments in the fu- 
ture to be able to check this prediction. In particular, 
we contrast it with a recent theory^^ in which there is 
a sufficiently large third neighbor exchange J3 and the 
monoclinic distortion alters the first and second near- 
est neighbor interactions, leading to the observed mag- 
netic order. In our model, the magnetic order is stabi- 
lized just by the anisotropy of J 2l with one ferromag- 
netic and one antiferromagnetic direction. (Refi£ii& also 
includes the anisotropic couplings but the third neighbor 
hopping makes it possible to stabilize the order with- 
out ferromagnetic directions.) Note, while doing the spin 
wave calculations, we expand about the equilibrium spin 
state, and hence the biquadratic interaction effectively 
leads to weak (w) and strong (s) nearest neighbor bonds, 
J\s Jiw Thus, our spin wave dispersion is repro- 
duced by a model involving only quadratic spin couplings, 
where JiS 1 -S 2 -X 1 (S 1 -S 2 ) 2 => (J 1 -2iifi<SrSa))Si-Sa 
Therefore the bonds between parallel spins are effectively 
weaker than those between antiparallcl spins, similar to 
Refi&i£. The spin wave spectrum is obtained using the 
Holstein-Primakoff expansion (see e.gJ^). Figure [5] com- 
pares the dispersions one expects in the two models, the 
first with a strong J3 and ours with a ferromagnetic di- 
agonal coupling. Note, the upper band of the dispersions 
curves in the opposite direction along the V-axis, because 
of the ferromagnetic coupling 
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FIG. 5: Comparison of the spin wave dispersion for order 
stabilized by J3 versus by a ferromagnetic diagonal interac- 
tion. The dispersions are plotted as a function of (kx, ky) — 
(k x — ky, k x + ky) in the upper pictures. The other graphs show 
the dispersion in the Y direction, a) The spin wave spectrum 
for J 2X = J 2 y = l,Ji = 1.47, K = 0.265, J 3 = 0.4, which 
have the same effective J\ w s as in Ref£. b) This shows the 
spectrum for J 2 x = —Jiy = l,Ji = 0.7, K = 0.265, J3 = 0. 
The sign of J 2 y is flipped to describe a ferromagnetic cou- 
pling, while J3 is set to zero. (Ji's value is decreased relative 
to a) to stabilize the order.) Note that the upper band changes 
from a hill shape to a valley. 



E. Doping Induced Incommensuration 

The properties of Fe±+ y Te have been experimentally 
investigated as y is varied. Experimentally, Fe\ +y Te is 
found to have incommensurate spin order when y is large 
enough, with an incommensurate wavevector that devi- 
ates from (7r/2,7r/2) linearly with doping^. One of the 
effects of the excess iron, which is believed to be in the 
Fe 2+ state, is to electron dope the system by 2y elec- 
trons. Here we consider how this may be explained as a 
result of the increasing electron density 6 = 2y in the 
Yz orbital, which makes nearest neighbor hopping more 
important. This will lead to incommensurate order that 
is related to the doping level. 

To see this we will consider two dimensional motion 
of the doped electrons, which can occur only if hopping 
amplitudes other than t 2 are taken into account. These 
types of motion are limited by the wide stripes of the 
commensurate (§, § ) spin pattern, increasing the kinetic 
energy, according to the uncertainty principle. Such ki- 
netic energy effect has a stronger dependence on the rel- 
ative angle between a pair of spins when the angle is 
close to 90° than the interactions considered above, so it 
is able to distort the commensurate collinear state with 
wave vector (|, — ? ) and <f> — tt/2, into a "double spiral", 
a slightly twisted version of the same state, of the form 
Eq. ([5]), and with <j> « | and k s» §, even when the 
doping is weak. 

To understand the kinetic energy effect, we assume a 
uniformly varying classical spin configuration and mini- 
mize the energy of a single Hund's coupled electron, Eq. 
@, hopping in this background. The resulting kinetic en- 



ergy (KE min ) times electron density, KE min S is added 
to the magnetic energy Q , and this total energy is min- 
imized to obtain the phase diagram in Fig. [5] A short- 
coming of this phase diagram is that the spiral phases 
may not be stable against phase separation (see below). 

We will now focus for simplicity on determining the 
kinetic energy due to hopping via t' 2 . We checked also 
that the nearest neighbor t\ of Ref. 14 leads to a very 
similar phase diagram. The kinetic energy is computed 
from the hopping Hamiltonian H t > 2 = — t' 2 J2i dyziri+x-— 
y)d,Yz{Yi) + h.c, which favors ferromagnetic alignment 
along the X-diagonal. Since the Hund's energy is larger 
than the i's, the natural basis for spin states on each site 
are the states quantized along the local spin orientation 
of the magnet (6*;). If the Hund's energy is assumed to be 
infinite, we can assume that the electron state is always 
aligned with the local spin. Such a state, with momentum 
p is represented by the following electron wavefunction, 

*=M e ;f ) e<p '"- <7) 

We can now find the variational energy of this wave func- 
tion, and minimize with respect to p. 

The expectation value of the kinetic energy for a bond 
between a pair of sites is 2t$l(ipl'ipj). Aside from the 
phase factor from the momentum, this overlap is pro- 
portional to cosfc for spins adjacent along the X direc- 
tion (the cosine of half the angle between the spins). 
In particular, because the spin is conserved during the 
hopping, the electron cannot hop at all onto a site 
with an antiparallcl spin. The electron kinetic energy is 
2t 2 cos(p x + p y ) + 2t' 2 cos( Px ~ Pv ) cos k. Minimizing over 
p, we obtain KE m i n = —2\t 2 \ — 2\t' 2 cos(k)\. Combin- 
ing this with the previously obtained energy in the ab- 
sence of t' 2 , Eqn. |6|), and minimizing with respect to 
<j) and k we obtain the double spiral. This can be eas- 
ily seen when the doping is very weak: focusing on the 
competing terms J 2 x and t' 2 5, and assuming k w | to 
find E = J 2X cos{2k) - 2i' 2 1 cos fc| » ~J 2X + 2J 2X (k - 

■|) 2 — 2t 2 5\k — ||, so \k — j\ = jj^ minimizes the en- 
ergy. The incommensuration therefore is proportional to 
the doping strength, as observed experimentally. 

If we allow the electrons to hop virtually into states 
that violate Hund's rule, the commensurate state sur- 
vives over a finite range of dopings. When the twist- 
ing is very small, it becomes easier for an electron that 
is trying to hop along the antifcrromagnetic diagonal to 
hop through the Hund's violating states than to stay in 
states that are parallel to the local spin but which have 
very small overlaps. In this case, there is no reason for 

the spin order to distort at all. The virtual hopping has 

t' 2 

the greater efficiency when \t' 2 cos k\ < |^-|, or (using the 
result for k at small 5), S < (A more detailed calcu- 
lation gives the same result.) Quantum fluctuations also 
allow the electrons to hop more easily in this direction, 
by temporarily making adjacent atoms become oriented 
parallel to one another. 
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Single Spiral 

FIG. 6: The phase diagram as a function of K and doping 
for Ji = \ JiXi t'2 = 1.5J2X, using the infinite Hund's cou- 
pling approximation. The three phases are the doubled spi- 
ral, the single spiral, and ordinary (ir, n) antiferromagnetism. 
Note that both the single and doubled spirals are constant 
along wavefronts parallel to Y. In the perpendicular direc- 
tion, the double spiral alternates between rotating through 
small angle and angles close to 180°, while the single spiral 
rotates through around 90° each time. A calculation which 
takes fluctuations or the finiteness of the Hund's constant into 
account would find a finite range of dopings where the com- 
mensurate phase is stable. The solid/dashed line represents a 
discontinuous/continuous transition. The discontinuous tran- 
sition boundary should actually be replaced by a wide swath 
of phase separation. 



A very similar phase diagram appears for t\ hopping. 
This hopping probably has weaker effects, however, since 
the orbital switches with each hop. This implies that 
an electron alternates between Hund's rule violating and 
satisfying states, so that the amplitude for a pair of these 
t 2 

steps is of order -j^. This process leads to an antiferro- 
magnetic interaction along rows and columns, but this 
causes incommensuration just as the ferromagnetism in- 
duced by motion along the diagonals did. 

Although these calculations explain how incommensu- 
rate order can occur, the order differs from the experi- 
mentally proposed pattern irA The proposed order has 
the spin dominantly along the Y direction whose magni- 
tude is modulated with an incommensurate wavevector, 
and also a spin spiral composed of the orthogonal spin 
directions, at the same wavevector. While further ex- 
perimental work is required to confirm the true nature 
of the complex incommensurate order, we note one in- 
teresting measurement, that orthorhombic symmetry is 
recovered at higher doping, (eg. at x — 0.1415,) where 
the incommensurate state is stabilized. This is consis- 
tent with our assumption that the monoclinic part of 
the distortion is strongly coupled to spin order. Fig. 
Q] shows that the sense of the monoclinic distortion is 
correlated with the bond energies Si ■ Sj. Shifting the 
spin ordering and thus the pattern of bond energies over 
one site would cause the lattice to tilt in the opposite 
direction along the Y axis. For an incommensurate or- 
der of the kind proposed in Ref4, these bond energies 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 
Doping 



FIG. 7: Energetics of phase separation. The solid curve shows 
the energy of states with uniform doping, with repulsion V = 
6J2X (upper curve) and without. The parameters are the 
same as in Fig. [6] with K — 0.4. The concave down portion 
of each energy curve is unstable to phase separation. The 
points labelled C show the unphysical transition between the 
doubled- and single-spiral states. The dashed line corrects 
the upper energy-curve, by showing the energy of the state 
where the A and B phases coexist. 

are also modulated with an incommensurate wavevector 
which removes the monoclinic distortion. In contrast, 
for the orders described by EqnfSl including both the 
double and single spirals we have considered, the bond 
energies are independent of incommensuration, hence a 
monoclinic distortion is expected throughout. An exam- 
ple of a spiral order that would not induce a monoclinic 
distortion, close to the commensurate state of interest is: 
Q i = {-l) x '+y>k{ Xl - yi), with k close to § . Effects that 
we have neglected, including spin anisotropy and the ef- 
fect of the excess iron moments, can modify the precise 
form of the incommensurate state. 

Phase Separation An additional shortcoming of this 
explanation of the incommensurate order is that the spi- 
ral phases can be unstable to phase separation^. The 
doped electrons prefer to be segregated into high and low 
density regions, with different spin orders. The predic- 
tions above can still be relevant, because repulsive forces 
between electrons help to limit the phase separation. 

Fig. [7] shows the energy as a function of doping both 
without and with a strong short-range repulsion, de- 
scribed by adding a term \V5 2 to the energy functional. 
In the figure, K = 0.4 and the other parameter values 
taken from Fig. [6] A concave down portion is unstable 
to phase separation. For a doping in between A and B, 
i.e., 0.1 < 5 < 0.6, dividing the compound up into two 
regions with the dopings 0.1 and 0.6 gives a state with a 
lower energy, whose energy is represented by the tangent. 
In this case the phase mixture is between a doubled spiral 
with a low electron density and an antiferromagnetically 
ordered high density portion. 
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The Coulomb interactions can have one of two conse- 
quences; they either stabilize a uniform state as just de- 
scribed (if the force is strong enough at short distances), 
or else they force the two coexisting phases to fill small 
alternating portions of the FeTe, rather than becoming 
completely segregated. The ordering with alternating re- 
gions of the (§,■§) and (tt, tt) phases could also be incom- 
mensurate. 

Additional Doping A surprising order may arise if the 
compound is sufficiently doped, which directly connects 
to our central assumption that the orthorhombic part of 
the structural distortion arises from orbital, rather than 
spin ordering. In our model, large electron doping would 
drive ferromagnetic order on each of the A and B sub- 
lattices. This optimizes the electron kinetic energy The 
sublattices continue to couple together antiferromagnet- 
ically, leading to (tt, tt) magnetic order. In Fig. [6l the 
order varies from commensurate, to the rhythmic spiral 
phase to the (tt, tt) antiferromagnetic phase as the dop- 
ing increases along a line at K = 0.5J2X- For a smaller 
value of K, the pure spiral phase occurs before (tt, tt) 
antifcrromagnetism sets in. Note however, that while 
(tt, tt) antiferromagnetic order is compatible with tetrag- 
onal symmetry, here we expect orbital order to persist 
and result in an orthorhombic distortion, clearly signal- 
ing the independence of structural and magnetic order. 

On the other hand, hole doping weakens the kinetic 
energy effect to the point where each sublattice has or- 
dinary antiferromagnetic order, and the biquadratic in- 
teraction orients the two sublattice-spins parallel to each 
other, leading to (tt, 0) magnetic order of the type seen 
in the FeAs materials. Again, while the (tt, 0) ordering 
would be accompanied by compression along the x or 
y axis, if this were magnetically driven, here we expect 
the orbital ordering to remain the same, so the diago- 
nal distortion would not change. A more direct test of 
whether the distortion or the magnetic ordering is the pri- 
mary phenomenon is to apply a strong magnetic field to 
eliminate the magnetic ordering and see whether the dis- 
tortion remains. This would be possible if the exchange 
interactions are weak, or can be weakened by modifying 
the compound somehow, e.g., by applying pressure. 



III. TRANSPORT AND SINGLE ELECTRON 
PROPERTIES 

The main consequence for transport of the orbital or- 
dering induced quasi one dimensionality is that the con- 
ductivity should be strongly anisotropic. Below the or- 
bital and magnetic ordering temperature, the excess elec- 
trons in the Y z band move much more readily in the 
Y direction, because is greater than ti,t' 2 . Further- 
more, the spins are ferromagnetically aligned along this 
direction, and can hence propagate easily. To travel in 
the orthogonal direction, they must cross through diag- 
onals where the spins are oppositely aligned. While a 
difference in electrical conductivity along the two diag- 



onal directions is hardly surprising given the symmetry 
of the low temperature state, the specific prediction here 
is that this will be a significant effect, and the nature of 
the anisotropy is that the low conductance is to be found 
along the antiferromagnetically ordered diagonals. 

Experimentally, a Drude peak has been observed to 
develop below the ordering transition 5 . Since a test of 
anisotropy demands a single domain, optical conductivity 
on a sample where the domain size is larger than the spot 
size is required, and should be feasible. There, depending 
on the direction of the polarization, a different conduc- 
tivity should result. Note, writing ct(uj) — ^^tt — , 
both the effective mass to* and the scattering rate r are 
expected to be anisotropic, which can be separately in- 
ferred from a knowledge of The anisotropy in this 
scenario is expected to be particularly prominent in the 
low frequency limit, since scattering by spins will limit 
the scattering time along the antiferromagnetic diagonal. 
At higher frequencies, only the anisotropy of to* will have 
an effect. However, if the scattering rate for the X direc- 
tion is extremely large, o~(w) will be flat in that direction, 
while conductivity along Y will have a Drude peak. 

ARPES Given the quasi one-dimensional dispersion 
along the diagonals, a narrow elliptical Fermi surface 
tilted at 45° to the Fe-Fe bonds, is expected to appear be- 
low the ordering temperature. Moreover, these would be 
orbitally polarized, which can be experimentally tested 
using polarized light to determine orbital content along 
high symmetry directions. When the scattering plane 
containing the photon and the ejected electron is per- 
pendicular to the surface, and along, say, the Xz plane, 
selection rules imply that S (P) polarized light with po- 
larization perpendicular to (parallel to) the scattering 
plane, ejects only electrons in the Yz (Xz) orbital. The 
strongly dispersing direction should therefore appear for 
S polarized light, and will have the form shown in Figure 
if a single domain is imaged. The polarization de- 
pendence provides an experimental signature even from 
multi-domain data. Current ARPES data on FeTe 4 has 
not reported such a signature, however, the intensity as- 
sociated with this Fermi surface segment is hard to es- 
timate reliably. Direct experimental probes of orbital 
ordering should also be able to test this scenario. We 
hope this prediction will stimulate further experiments. 



IV. CONCLUSIONS 

We have argued that FeTe, a material closely related 
to the recently discovered Fe based superconductors, is 
likely to be a fairly strongly correlated material. This mo- 
tivates us to use a local picture of the electronic structure, 
which in turn can explain, quite naturally, the unusual 
magnetism observed in this material. The key ingredient 
of this scenario is an emergent quasi one dimensionality 
arising from orbital ordering along the diagonals. While 
this theoretical scenario is a simplified caricature of the 
real system, it does make several qualitative predictions 
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for experiments which should be readily testable. The 
conductivity should be anisotropic, with higher conduc- 
tivity along the ferromagnetic direction, if the unusual 
magnetic order is caused by the kinetic energy of con- 
duction electrons. Certain parts of the electron Fermi 
surface should disappear in polarized ARPES, indicat- 
— ing the orbital ordering. Lastly, if the Jahn- Teller effect 

k x causes the orthorhombic distortion, then the distortion 

should persist even when the magnetic order is removed 
or changed to one which would not be expected to favor 
lattice distortion on the basis of symmetry. 



FIG. 8: The shape of the Fermi surface of the Yz electrons 
in the actual Brillouin zone. We predict that the surface is 
narrow in the ferromagnetic direction. Furthermore, certain 
parts of the ellipse should disappear in ARPES with polar- 
ized light. The dashed portions of the Fermi surface fade 
out for S polarization. This tests the assumption that the 
anisotropy results not from the lattice distortion but because 
the dispersing electrons are in a single orbital. 
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